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Abstract. We analyze the problem of parameter estimation for compact binary 
systems that could be detected by ground-based gravitational wave detectors. So 
far this problem has only been dealt with for the inspiral and the ringdown phases 
separately. In this paper, we combine the information from both signals, and we 
study the improvement in parameter estimation, at a fixed signal-to-noise ratio, 
by including the ringdown signal without making any assumption on the merger 
phase. The study is performed for both initial and advanced LIGO and VIRGO 
detectors. 



PACS numbers: 04.80.Nn, 95.55.Ym, 97.60.Gb, 07.05.Kf 
1. Introduction 

Coalescing compact binaries consisting of either black holes (BH) or neutron stars 
(NS) are among the targets of on-going searches for gravitational waves in the 
data of ground-based interferometric detectors such as GEO 600 |2] , the Laser 
Interferometer Gravitational Wave Observatory (LIGO) [3 i], TAMA 300 and 
VIRGO 0. 

The coalescence of a compact binary system is commonly divided into three 
stages, which are not very well delimited one from another, namely the inspiral, the 
merger and the ringdown. Many studies so far have focused on the gravitational waves 
emitted during the inspiral phase because the inspiral waveform is very well understood 
13 m 13 Uni El 113 nil IH and the event rates seem promising [H UU El HHI UH . 
Gravitational waves from the merger can only be calculated using the full Einstein 
equations. Because of the extreme strong field nature of this epoch neither a 
straightforward application of post-Newtonian theory nor perturbation theory is very 
useful. Recent numerical work [201 1^11 122L 1231 has given some insights into the 
merger problem, but there are no reliable models for the waveform of the merger phase 
at this time. The gravitational radiation from the ringdown phase is also well known 
and it can be described by quasi-normal modes |25) . In spite of the importance of the 
ringdown there are a fewer publications on ringdown searches compared to those for 
inspiral searches. 

Flanagan and Hughes were the first in studying the contribution of the 
three phases to the signal-to-noise ratios both for ground-based and also space-based 
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interferometers, but they did not study the problem of accuracy in the parameter 
estimation. This is an important problem because many efforts are now underway to 
detect both the inspiral and the ringdown signals using matched filtering techniques 
in real data [271 128L 1291 13(J[ 13 1| . In a recent paper parameter estimation of 
inspiralling compact binaries has been revised using up to the 3.5 restricted post- 
Newtonian approximation, extending previous analysis |33l I34j . but ignoring the other 
stages. The parameter estimation for the ringdown phase alone has also been studied, 
some time ago, for ground based detectors [SH1ES1E21 as well as for LISA [SH|- The 
aim of this paper is to discuss how parameter estimation can be improved by using 
information from both the inspiral and the ringdown phases combined together in 
matched filtering like analysis for different ground-based detectors. 

This paper is organized as follows: Section |21 introduces our notation and reviews 
the basic concepts of signal parameter estimation in matched filtering. Section 13 
provides the noise curves used in this study for initial and advanced LIGO and 
VIRGO. Section ^ briefly describes the waveforms that we are looking for. For the 
inspiral phase, we consider a non-spinning compact system with circular orbits and 
the waveform in the restricted post-Newtonian approximation. For the ringdown, we 
assume that the dominant mode has I = m = 2 and therefore the waveform is given 
by an exponentially decaying sinusoid. Section [3 studies the impact on the parameter 
estimation for coalescing binary black holes, by combining the signals from both the 
inspiral and ringdown phases and compares the results with the case of inspiral phase 
alone. The results are presented for a fixed inspiral signal-to-noise ratio of 10. Different 
number of parameters are used as well as different values for the ringdown efficiency. 
Finally section concludes with a summary of our results and plans for further work. 
In the Appendices we collect various technical calculations and we present an explicit 
analytical calculation of the Fisher matrix for the ringdown phase that has been used 
to compare with the numerical results. 



2. Summary of parameter estimation 

In this section we briefly review the basic concepts and formulas of signal parameter 
estimation relevant to the goal of this paper; we refer the reader to |33| for a more 
detailed analysis. 

The output of a gravitational wave detector can be schematically represented as 

s{t) = h{t) + n{t) , (1) 

where n{t) is the noise that affects the observation and h{t) is the gravitational 
wave signal measured at the detector, a linear superposition of the two independent 
polarizations of the strain amplitude /i+ and hx, given by 

h{t) ^ F+{e,cj),^P)h+{t) + FxiO,(^,ij)hx{t) , (2) 

where and -Fx sltc the antenna pattern functions, that depend on the direction 
of the source in the sky {9, (f) and the polarization angle V'- In case of a laser 
interferometer detector, the expressions of and F^ are given by |44j : 

F+{0,(f>,ij) = i(l + cos2e')cos20cos2-i/'-cos6'sin2(/«sin2V' , (3) 
Fx{0,(j),ij) = i(l + cos^6')cos20sin2?/>-|-cos6'sin20cos2V' . (4) 
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For sake of simplicity wc shall made the standard assumptions that the noise 
n{t) has zero mean and it is stationary and Gaussian, although in realistic cases this 
hypothesis is likely to be violated at some level. Within this approximation, the 
Fourier components of the noise are statistically described by: 

E[hif)h*U')]^hif-r)SM), (5) 

where E[] denotes the expectation value with respect to an ensemble of noise 
realization, the * superscript denotes complex conjugate, 5'„(/) is the one sided 
noise power spectral density, and tildes denote Fourier transforms according to the 
convention 

/oc 
e*2-/*x(t) dt . (6) 
-oo 

With a given noise spectral density for the detector, one defines the "inner 
product" between any two signals g{t) and h{t) by: 

With this definition, the probability of the noise to have a realization rig is just: 

p{n = no) oc e(""l"")/2 . (8) 

The optimal signal-to-noise ratio (SNR) p, achievable with linear methods {e.g., 
matched filtering the data) is given by the standard expression 

In the limit of large SNR, if the noise is stationary and Gaussian, the probability 
that the gravitational wave signal h(t) is characterized by a given set of values of the 
source parameters A = {A'^jfe is given by a Gaussian probability of the form |35| : 

1 

where AA*"' is the difference between the true value of the parameter and the best- 
fit parameter in the presence of some realization of the noise, p^"-* (A) represents the 
distribution of prior information (a normalization constant) and Tjk is the so-called 
Fisher information matrix defined by 

Tij = {^^h\^Jh) = 2 -^-^^ df , (11) 

where di = 

The inverse of the Fisher matrix, known as the variance-covariance matrix, gives 
us the accuracy with which we expect to measure the parameters A'^ 

^ (Y-iyk ^ (AA^AA'^) . (12) 

Here the angle brackets denote an average over the probability distribution function 
in Eq. (|10|l . The root-mean-square error ak in the estimation of the parameters A'^ 
can then be calculated, in the limit of large SNR, by taking the square root of the 
diagonal elements of the variance-covariance matrix, 

au^{(^\^fY'^^VYF, (13) 
and the correlation coefficients c^^ between two parameters A^ and A*^ are given by: 



p{\\h) =p(°)(A)exp 



-;r,,AA^AA^ 



(10) 
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3. Noise spectra of the interferometers 



In this paper, we use three different noise curves to understand the effect of detector 
characteristics on the parameter estimation. The noise curves used are initial and 
advanced LIGO and VIRGO as in Those are: 
For the initial LIGO 

So [(4.49a;)-5S + 0.16x-'^-^^ + 0.52 + O.S2x^] , / > / 
oo , f < 



Snif) 



(15) 



where 



f/fo, with /o = 150 Hz (a scaling frequency chosen for convenience), 



fs = 40 Hz is the lower cutoff frequency, and S'o = 9 x 10 "'^ Hz 
For advanced LIGO the noise curve is given by 



^n(/) 



So 



5x 



lll(l-.T^+x''/2) 



OO 



f>fs 
f<fs 

where fo = 215 Hz, f, = 10 Hz and S'o = 10'^^ Hz'^. 

Finally, for the VIRGO detector the expected noise curve is given by: 



(16) 



Snif) 

where fo = 500 Hz, fs 



20 Hz and So = 3.24 x lO"''^ Hz"!. 



So [(6.35a;)- 
oo , 



f>fs 
f<fs 



(17) 



4. The gravitational-wave signal 

As discussed in the introduction, the coalescence and its associate gravitational wave 
signal can be divided into three successive epochs in the time domain: inspiral, merger 
and ringdown. During the inspiral the distance between the stars diminishes and the 
orbital frequency sweeps up. For low-mass binary systems, the waveforms are well 
modeled using the post-Newtonian approximation to general relativity |71 El . 
Eventually the post-Newtonian description of the orbit breaks down, and the black 
holes cannot be treated as point particles any more. What is more, it is expected 
that they will reach the innermost stable circular orbit (ISCO), at which the gradual 
inspiral ends and the black holes begin to plunge together to form a single black 
hole. This is referred as the merger phase. At present, the merger phase is not well 
understood and no analytical reliable waveforms exist. At the end, the final black 
hole will gradually settle down into a Kerr black hole. The last gravitational waves 
will be dominated by the quasi-normal ringing modes of the black hole (see |41| and 
references therein) and can be treated using perturbation theory 0^1 • At late time, 
the radiation will be dominated by the I ^ m = 2 mode [231 • This is the so-called 
ringdown phase. 

The gravitational waveform of coalescing compact binaries thus takes the form 

{/li„=pi™i(i) -oo < t < Tisco 
/l„o,,o,(i) Tisco < < < T^^^ (18) 
/lri„Edow„(i) Tqnr < i < oo 

where Tjsco is the time when the system reaches the ISCO and Tq^r is the time when 
the quasi-normal mode I = m = 2 begins to dominate the ringdown, although there 
is some arbitrariness in choosing Tjsco and Tq^r to delimit the three epochs. 



Parameter estimation of compact binaries 



5 



The inspiral waveform 

For a non-spinning compact binary system with circular orbits, the two polarizations 
/i+ and hx of the inspiral waveform can be well described by the post-Newtonian 
expansion. Thus setting G = c ~ 1, they read: 

, (19) 

where v = {Muj)'^/'^, lo is the orbital frequency, r is the distance to the source, 
M = mi -|- 7712 is the total mass, /i ~ mim2/M is the reduced mass, rj — fi/M is 
the symmetric mass ratio and A4 = /i^/^M^/^ = Mrf^^ is the chirp mass. 

In what follows we consider the waveform in the restricted post-Newtonian 
approximation |48| . corresponding to a frequency twice the orbital frequency, and 
we ignore higher order harmonics. This corresponds to the lowest terms in the series 
(HHl. The functions iJ^"' are given by [101 : 

= - (1 + cost2)cos$(t) , (20) 

iJ^"^ = - 2 cos i sin$(t) , (21) 

with b being the angle between the orbital angular momentum of the binary and the 
line of sight from the detector to the source. $ is the phase of the gravitational wave an 
the instant t, that we consider modeled through 2nd post-Newtonian order, neglecting 
the higher order terms in this analysis, since they would not contribute significantly 
to the result. 

The Fourier transform of the inspiral waveform can be computed using the 
stationary phase approximation [2311^123001 and this yields: 

r (r,l A.sr^/V*(/) /</:sco 

with 



Ans = -^^J^^-"'^Fl(l^^^Y^Fl^\ (23) 

*(/) = 2^/<, - 0, _ J + _ ^ Aku''-^ , (24) 

where c — cos i, tc refers to the coalescence time, 0c is the phase at the coalescence 
instant, u — (ttA^/)^/'^, and the coefficients Ak are given by 

^0 = 1 (25) 

^1=0 (26) 

20,743 11 , _2/5 

M = — ( \ -11)11 ' (27) 

9 ^336 4 ' ^ ' 

A3 = - 167r?r^/^ (28) 

,3058673 5429 617 , , 

= 10^1^6064 + 1008^ + T44'^>^' (^^^ 

We also consider the ISCO to take place at a separation of 6Af , corresponding to a 

final frequency 

= ■ ^^^^ 
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4-2. The ringdown waveform 



The ringdown portion of the gravitational wave signal we consider can be described as 
the I = m ~ 2 quasi-normal mode. Therefore the gravitational radiation in the time 
domain is expected as the superposition of two different damped sinusoids, although 
one of these exponentials could be invisible in the actual waveform as discussed in |88j . 
In our study, we assume that the ringdown waveform can be written as in j25j . 
corresponding to a circularly polarized wave. In this way we have 



h+{t)-ih^{t) = —^Si{i,f3,a) 



X exp 



-i27r/QNR(i - to) 



'-{t - to) + m) 



(31) 



where to is the start time of the ringdown, ipo the initial phase, M is the total mass of 
the system mass (see 0^1 for further discussions), /qnr and Q arc the central frequency 
and the quality factor of the ringing. For this mode, a good fit to the frequency /qnr 
and quality factor Q, within an accuracy of 5%, is 

/qnr 0.63(1 -a)3A"]^, (32) 

Q «2(l-a)-9/20^ (33) 

where aM^ is its spin, and a is the Kerr parameter that lies in the range (0.0,0.998). 
In our study we set a to the near extremal value of 0.98 (as in I^HDj although we 
consider a as any other independent parameter when evaluating the Fisher matrix. 
The function 2S2 is the spin weighted spheroidal harmonic that depends on the 
inclination angle of the black hole axis seen from the observer and the Kerr parameter 
a. A is a dimensionless coefficient describing the magnitude of the perturbation when 
the ringdown begins. Although the value of the amplitude is uncertain, we set the 
amplitude of this mode by assuming that a fraction e of the system's mass is converted 
into gravitational waves during the ringdown |26| 

2 



Therefore 



A 



1 



-A^mVqnrQ = eM 



128 7f e 



, Mf^^nQ 

The strain produced at the detector can be written as 

{t - to) ''■/qnr 



(34) 



(35) 



/iRD(t) = Ard exp 



Q 



where 



AM 



The Fourier transform of the waveform becomes: 

ZTT 



cos(-27r/QNR(t - to) + 70) (36) 

(37) 
(38) 



^'q'" ~ 2i(/ — /qnr) '''q'" ~ 2i(/ + /qnr)^ 
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5. Parameter estimation of compact binaries using the inspiral and 
ringdown waveforms 

In this paper we want to study the impact on the parameter extraction by combining 
the signals from the inspiral and the ringdown epochs, neglecting all information 
coming from the merger epoch itself since no reliable waveforms exist so far. 

Following earlier works, we choose the set of independent parameter Aj^g 
describing the inspiral signal to be 

Aims = {In Ans, /o^c, 0c, InTW, Inr?} , (39) 

while for the ringdown, the parameters could be 

Ard = {lnylR.D,lnM,lna,7o,to} . (40) 

A possible approach to this problem of parameter extraction would be to consider 
the two sets of parameters ()39|l and (|40|l as independent, perform matched filtering 
using the two template families (for the inspiral and the ringdown waveforms) and then 
reduce the uncertainties in the parameter estimation (in particular for the masses) by 
making a posterior consistency check |46) . However in this paper we follow a different 
approach. Wc consider only a single coalescing waveform, as if we were performing 
matched filtering with a single template family bank, as given by Eq. H18() . which 
describes the different phases, but ignoring the information from the merger phase, 
and we focus our attention on how parameter estimation of the Ams parameters can 
be improved as compared to the case in which the inspiral waveform is used alone. For 
this reason the study presented here focuses only in those mass ranges for which the 
inspiral signal alone could be detectable by the detector. This corresponds to a total 
mass of approximately 1-100 Mq for initial LIGO, 1-400 Mq for advanced LIGO, 
and 1-200 Mq for VIRGO. 

The global waveform considered here becomes in the Fourier domain 

hodf) = k.s{.f) + Kn{f), (41) 

where ft.iNs(/) and /ird(/) are given by equations and respectively. This 
global waveform is completely determined by a set of independent parameters, given 

by 

Agl = {lnANs,/oic,<?!'c,ln7W,ln77,lna,7o,io} • (42) 

Notice that we do not include Iu^Ird as an independent parameter since for a given 
source location and orientation, and a given ringdown efficiency e, the ringdown 
amplitude A RD 5 given by Eq. (j37|l . is determined by the inspiral amplitude yli^s, 
the Kerr parameter and the masses. Instead, what we do is to find a heuristic relation 
between A^tu and Ams by averaging over source directions and black-hole orientations, 
making use of the angle averages: = {Fx)e,<p,i, = 1/5, (F+Fx)e,0,^ = 0, 



1/3, ((1 + c2)2), = 28/15, and (|2.S|P).,/3 = l/47r. The angle averaged root 
mean square (rms) values of the inspiral and ringdown amplitudes Amsi A-ru, given 
by Eq. ^ and become 

A^l ^ = -7^^- (43) 
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From the above equations, together with Eq. Ij^Sfl . we derive a relation A^o — 
^rd(^ins5-^i tt) through the ratio of rms of both amphtudes: 

Ao(Ans, X, 7/, a) = 4^ Ans = \ TrT^'^ 'l'^'^'^' . (45) 

A'jTs V ^^/Qnh V 

Note that in this relation the product of M and /qnr is just a function of a as can be 
seen from Ea. (|32|l . This relation (|45|1 will be used in calculating the SNR as well as 
the Fisher matrix for the global waveform. 

The SNR values for equal-mass black hole binaries are shown in figure ^ This 
figure indicates that there is a range of masses (different for the different noise curves) 
for which both the inspiral and the ringdown signals could be detectable and one could 
search for both portions of the signal in order to improve the SNR and the accuracy 
in parameter estimation^. 

It is clear that both in the time domain, as well as, in the frequency domain, the 
inspiral signal is decoupled from the ringdown one. The inspiral waveform /iins(/) 
ranges from the lower cut-off frequency fg to /iscoj while the ringdown /Ird(/) is 
centered around /qnr with a certain bandwidth, that in the literature is considered 
to be smaller than A///qnr = 0.5§, although in our numerical simulations for the 
ringdown signal we use the same lower cut-off frequency fs and a higher cut-off 
frequency of 5000 Hz. This justifies that the Fisher matrix, defined in Eq. Hll|) . 
of the global waveform can be computed as 

Tij ^ (3i/iGL|9j7iGL) = idihms\djhms) + {,dih^,o\djh^D) , (46) 

neglecting the cross elements (9i/iiNs|3jft.RD)- Therefore the Fisher matrix can be 
computed as the sum of Fisher matrix of the inspiral waveform plus the Fisher matrix 
of the ringdown 

Tgl = r™s + , (47) 

where we just need to be consistent in computing the elements corresponding to the 
same parameter set. Also the total SNR is given by 

Pgl = Pms + Prd ■ (48) 
The way we proceed is to analyze first the well known case of the inspiral 
signal alone, and then we compare the results with those when using the inspiral 
and ringdown waveforms together. In order to separate the effects of increasing the 
number of parameters from the fact we are using a more complex waveform, we study 
two different cases: 

i. The case in which only the five inspiral parameters are considered. This is 
equivalent to have no uncertainties in the spin of the final black hole, nor in 
the initial phase and time of the ringdown signal. This of course, would not 
be realistic in a search, but it provides the optimal improvement in parameter 
estimation one could expect from the fact that we added the ringdown waveform. 

X The calculation of the SNR for the ringdown waveform is computed differently from what was done 
by Flanagan and Hughes in I2bl . Instead of taking |t — to| in the damped exponential, integrating 
over t over — oo to -|-oo and dividing the result by \/2 to compensate for the doubling, we assume 
that the waveform /iRD(t) vanishes for t < tg and integrate only over t > tg. 

§ Note that the distance between /qnr and /isco is larger than the bandwidth of the ringdown 
signal /ird(/). In particular if we consider the value a = 0.98 then (/qnr — /isco)//qnr = 0.833 
that suggests no overlap between the inspiral and the ringdown signal. 
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Figure 1. The averaged signal-to-noise ratio for equal-mass black hole 
coalescences detected by ground-based interferometers at a luminosity distance 
of 1 Gpc. The solid line is the SNR curve for the inspiral, and the dash and dash- 
dotted lines for the ringdown portion of the signal assuming a value of a = 0.98 
and e equals to 1.5% and 0.5%, respectively. The top panel corresponds to initial 
LIGO, the middle panel to advanced LIGO and the bottom one to VIRGO. 
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ii. The more realistic case in which ah the eight independent parameters given in ()42|) 
are considered. 

In [Appendix A| the reader can find the exphcit calculations of all the waveform 
derivatives necessary to compute the Fisher matrix which is then computed 
numerically. 

In our analysis, we set the Kerr parameter a — 0.98 (as in HHDj and we consider 
two different values of e: a more optimistic value of 1.5%, (half the value used in pHjl. 
and a more pessimistic one of 0.5%, that are more consistent with recent numerical 
simulations With these parameters we study a range of masses, analyzing both 
the equal-mass and unequal-mass cases, for three ground-based detectors: initial 
LIGO, advanced LIGO and VIRGO, using the noise curves described in section |21 
All the errors are computed at a fixed value of inspiral SNR of 10. 

For the equal-mass case the results are presented in figures |21 El and ^ 
corresponding to initial LIGO, advanced LIGO and VIRGO, respectively. The errors 
of tc, -M and rj and some of the associated correlation coefficients for the inspiral signal 
alone, for different pairs of masses, can be found in table ^ In tables El and |21 one 
can find the comparison of errors and correlation coefficients for the different cases we 
have analyzed. In particular table|21refers to the case (i) in which only the five inspiral 
parameters are considered and table 13 refers to the case (ii) in which we use the eight 
global parameters. In all cases the errors improve, and the improvement is higher for 
larger masses for which the ringdown signal contribution to the SNR increases. This 
improvement could be explained by the greater structure and variety of the global 
waveform but also by the variation of some of the correlation coefficients, although 
this is not fully assessed in this paper. We have just noticed that the correlations 
coefficients relative to the masses decrease when the ringdown signal is added, as can 
be seen in the tables. We also notice that the improvement is very significant for 
massive systems with very large errors for the inspiral waveform alone. These large 
errors are associated to the small number of useful cycles of the inspiral signal of 
these systems [32]. Therefore the effect induced in parameter estimation due to the 
inclusion of the ringdown signal could be understood in terms of additional number of 
gravitational wave cycles accumulated. Although, from the present analysis, it is not 
clear which of these considerations is the dominant aspect to completely understand 
the variation in parameter estimation observed with the global waveform. 

The numerical results for the inspiral and the ringdown waveforms separately 
have been verified by comparing with those existing in the literature (for different 
masses and noise curves). Moreover, for the ringdown case alone we have also found 
a good agreement with an analytical approximation as described in [Appendix B| 

6. Conclusions 

We have carried out a study to understand the implications of adding the ringdown 
to the inspiral signal on parameter estimation of non-spinning binaries using the 
covariance matrix. We have compared the results using three different noise curves 
corresponding to initial LIGO, advanced LIGO and VIRGO. 

The result of our study is that the parameter estimation of tc, 4>c, M. and rj 
improves significantly, as expected, by employing the extra information that comes 
from the ringdown, for those systems with a total mass such that both the inspiral and 
the ringdown signal could be detectable by the detectors. Naturally the improvement 
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Figure 2. In this figure we compare the errors in the estimation of ic, M and 
rj for equal-mass black hole coalescences by the initial LIGO interferometers at a 
fixed inspiral SNR of 10. The solid line corresponds to the inspiral signal only and 
the others to the combined inspiral plus ringdown waveforms. The dashed line 
corresponds to the case in which only the five independent inspiral parameters 
I39i are considered for e = 1.5%, while the dot-dashed lines correspond to the 
cases in which we consider all the independent global parameters I42i and e equals 
to 1.5% and 0.5%, respectively. 
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Figure 3. Same as figure 1^ for advanced LIGO. 




Figure 4. Same as figure for VIRGO 
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Table 1. Measurements of errors and some of the associated correlation 
coefficients using the 2nd post-Nowtonian binary inspiral waveform at a SNR of 
10. For each of the three detector noise curves the table presents Ate (in msec), 
A(pc (in radians), AAi/A4 and Arj/r] (in percentages). The cases considered here 
correspond to NS-BH and BH-BH binaries of different masses. 



mi 


m2 


Ate 




AM/M 


Arj/rj 








(Mo) 


(Mo) 


(msec) 


(rad) 


(%) 


(%) 








Initial LIGO 


20 


1.4 


3.3708 


7.8897 


0.7873 


10.6451 


0.9294 


0.9760 


0.9847 


50 


1.4 


47.2786 


65.2176 


6.5073 


50.6737 


0.9845 


0.9953 


0.9967 


20 


10 


7.6032 


15.0441 


8.0579 


72.3184 


0.9573 


0.9851 


0.9918 


60 


10 


253.505 


291.931 


176.566 


1058.45 


0.9947 


0.9983 


0.9990 


Advanced LIGO 


25 


1.4 


1.5726 


2.3121 


0.0769 


2.0304 


0.8002 


0.9397 


0.9319 


100 


1.4 


20.5373 


11.5863 


0.3354 


0.9734 


0.9426 


0.9829 


0.9859 


200 


1.4 


171.714 


59.9235 


1.5215 


12.2235 


0.9836 


0.9951 


0.9965 


100 


10 


24.9199 


13.3081 


2.5514 


26.7760 


0.9500 


0.9843 


0.9886 


200 


10 


205.361 


69.9572 


12.1978 


92.3077 


0.9855 


0.9955 


0.9970 


100 


50 


58.8914 


25.7410 


17.1443 


133.132 


0.9707 


0.9901 


0.9943 


175 


50 


282.241 


93.6166 


66.6933 


434.669 


0.9884 


0.9962 


0.9978 


VIRGO 


20 


1.4 


1.5875 


2.8685 


0.1339 


3.0661 


0.8619 


0.9499 


0.9679 


100 


1.4 


89.3534 


62.0062 


3.1241 


24.8609 


0.9839 


0.9951 


0.9966 


20 


10 


2.9432 


4.3971 


1.1712 


16.7567 


0.9025 


0.9630 


0.9805 


50 


10 


14.9211 


14.5300 


4.7808 


46.7855 


0.9564 


0.9853 


0.9912 


100 


10 


128.487 


85.3035 


28.3745 


202.072 


0.9874 


0.9960 


0.9975 


70 


50 


198.798 


128.504 


135.517 


804.327 


0.9907 


0.9968 


0.9983 


90 


50 


509.639 


299.335 


338.459 


1865.97 


0.9948 


0.9982 


0.9991 



is larger in the case of considering a smaller number of parameters, but in both cases, 
the five parameter case and the eight parameter case, the improvement is significant. 
The study is performed at a fixed inspiral SNR of 10, therefore the error in A^ins/-4ins 
would be of 10% for the inspiral signal alone. This is also improved by adding the 
ringdown. 

In this work we have made a number of simplifying assumptions, ignoring the 
merger phase, considering only the 2nd post-Newtonian inspiral phase formula instead 
of the 3.5 that is already known, using only a single mode for the ringdown signal, 
and ignoring angular dependencies (because of the angle averages we use). For this 
reason the results obtained here should be considered just as an indication of which 
could be the real effect in the parameter estimation by combining the inspiral with the 
ringdown signal. The preliminary results obtained here seem to be very encouraging. 
Therefore it would be interesting to extend the analysis to a more realistic case and also 
for different data analysis techniques (different from matched filtering). Although we 
have focused on ground-based detectors a similar study could be performed for LISA. 
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Table 2. Measurements of errors and associated correlation coefficients using the 
2nd post-Newtonian binary inspiral waveform at a SNR of 10 together with the 
ringdown waveform, using a set of five parameters {In^ms, /o4c, 0ci InA^, In?;}, 
excluding {Ina, 7o,to}- 



e 


mi 


m2 


-4lNS 




A(j)c 


AM/M 


Av/V 


Ct^M 






Pgl 


(%) 


(Mo) 


(Mq) 


(%) 


(msec) 


(rad) 


(%) 


(%) 










Initial LIGO 


1.5 


20 


1.4 


9.9990 


3.3424 


7.8207 


0.7807 


10.5517 


0.9282 


0.9756 


0.9844 


10.0010 


1.5 


50 


1.4 


9.9697 


14.8477 


19.9568 


2.0873 


15.5191 


0.8468 


0.9513 


0.9680 


10.0304 


1.5 


20 


10 


9.9852 


2.5040 


4.5312 


2.6880 


21.6438 


0.6113 


0.8537 


0.9239 


10.0148 


1.5 


60 


10 


9.2208 


9.7683 


1.7464 


10.3674 


17.5337 


-0.9594 


-0.9369 


0.9922 


10.9880 


0.5 


20 


1.4 


9.9997 


3.3613 


7.8665 


0.7851 


10.6137 


0.9290 


0.9759 


0.9846 


10.0003 


0.5 


50 


1.4 


9.9899 


23.1830 


31.6986 


3.2146 


24.6364 


0.9362 


0.9803 


0.9866 


10.0101 


0.5 


20 


10 


9.9951 


3.7473 


7.1810 


3.9909 


34.4457 


0.8250 


0.9372 


0.9660 


10.0050 


0.5 


60 


10 


9.6889 


10.2048 


2.1410 


10.8297 


18.7519 


-0.9490 


-0.8894 


0.9796 


10.3398 


Advanced LIGO 


1.5 


25 


1.4 


9.9994 


1.5723 


2.3116 


0.0769 


2.0300 


0.8002 


0.9396 


0.9318 


10.0006 


1.5 


100 


1.4 


9.5820 


14.3057 


7.9691 


0.2351 


2.7291 


0.8825 


0.9644 


0.9711 


10.4363 


1.5 


200 


1.4 


7.9748 


27.1408 


8.3128 


0.2626 


1.7010 


0.4037 


0.7836 


0.8779 


12.5470 


1.5 


100 


10 


7.4090 


3.9747 


0.9300 


0.4705 


1.4721 


-0.6895 


-0.1562 


0.7144 


13.5008 


1.5 


200 


10 


4.4008 


15.1787 


1.2639 


1.2077 


2.1026 


-0.9234 


-0.8621 


0.9726 


23.1557 


1.5 


100 


50 


3.9049 


6.6557 


0.8229 


2.2776 


3.8235 


-0.8769 


-0.8672 


0.9951 


28.3348 


1.5 


175 


50 


4.7381 


16.8821 


1.1369 


5.3223 


8.8739 


-0.9333 


-0.9325 


0.9997 


47.1507 


0.5 


25 


1.4 


9.9998 


1.5725 


2.3119 


0.0769 


2.0302 


0.8002 


0.9397 


0.9319 


10.0002 


0.5 


100 


1.4 


9.8546 


17.6073 


9.8910 


0.2882 


3.3904 


0.9221 


0.9766 


0.9809 


10.1475 


0.5 


200 


1.4 


9.1638 


42.2225 


14.0443 


0.3882 


2.8678 


0.7397 


0.9155 


0.9452 


10.9152 


0.5 


100 


10 


8.8592 


4.2986 


1.2823 


0.4993 


2.3083 


-0.4740 


0.2010 


0.6867 


11.2882 


0.5 


200 


10 


6.4175 


15.3235 


1.4330 


1.2144 


2.2822 


-0.8994 


-0.7383 


0.9295 


15.6651 


0.5 


100 


50 


5.6355 


6.6795 


0.8354 


2.2862 


3.8925 


-0.8744 


-0.8461 


0.9858 


18.2835 


0.5 


175 


50 


5.2676 


17.0019 


1.1474 


5.3641 


8.9507 


-0.9338 


-0.9315 


0.9991 


28.4206 


VIRGO 


1.5 


20 


1.4 


9.9884 


1.5739 


2.8419 


0.1328 


3.0372 


0.8595 


0.9490 


0.9673 


10.0116 


1.5 


100 


1.4 


9.8123 


23.6168 


15.7433 


0.8518 


6.3184 


0.7771 


0.9281 


0.9533 


10.1916 


1.5 


20 


10 


9.8860 


1.3595 


1.8045 


0.5264 


6.7294 


0.5307 


0.8124 


0.8998 


10.1154 


1.5 


50 


10 


9.5906 


2.5435 


1.4330 


0.8642 


4.2425 


-0.4119 


0.2591 


0.7145 


10.4273 


1.5 


100 


10 


8.7416 


8.8760 


1.5912 


2.6407 


4.9586 


-0.8994 


-0.7351 


0.9318 


11.4628 


1.5 


70 


50 


8.0221 


10.4639 


1.1242 


9.8663 


16.4643 


-0.9298 


-0.9263 


0.9987 


15.1284 


1.5 


90 


50 


8.8177 


12.3529 


1.0874 


11.8368 


19.7325 


-0.9136 


-0.9119 


0.9994 


18.3197 


0.5 


20 


1.4 


9.9961 


1.5830 


2.8595 


0.1335 


3.0564 


0.8611 


0.9496 


0.9677 


10.0039 


0.5 


100 


1.4 


9.9362 


37.4447 


25.6283 


1.3239 


10.2789 


0.9094 


0.9719 


0.9808 


10.0643 


0.5 


20 


10 


9.9616 


1.8724 


2.6772 


0.7371 


10.1277 


0.7565 


0.9058 


0.9501 


10.0386 


0.5 


50 


10 


9.8577 


3.1084 


2.2677 


1.0361 


7.0748 


0.0358 


0.5871 


0.7987 


10.1444 


0.5 


100 


10 


9.5177 


9.1524 


2.1289 


2.6903 


5.9890 


-0.8220 


-0.4653 


0.8587 


10.5102 


0.5 


70 


50 


8.7160 


10.3890 


1.2208 


9.7697 


16.3850 


-0.9375 


-0.9251 


0.9955 


11.9564 


0.5 


90 


50 


9.0164 


14.0215 


1.1789 


13.7140 


22.8751 


-0.9326 


-0.9288 


0.9987 


13.3618 
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Table 3. Measurements of errors and associated correlation coefficients 
using the 2nd post-Newtonian binary inspiral waveform at a SNR of 10 
together with the ringdown waveform, using a set of eight parameters 
{In Ans, /otc, <f>c, In A^, In 77, In a, 70, to}- 



e 


mi 


7712 




Ate 


A(j,c 


AM/M 


Ar;/?7 










(%) 


(Mo) 


(Mo) 


(%) 


(msec) 


(rad) 


(%) 


(%) 










Initial LIGO 


1.5 


20 


1.4 


9.9991 


3.3699 


7.8874 


0.7871 


10.6420 


0.9293 


0.9760 


0.9847 


10.0010 


1.5 


50 


1.4 


9.9726 


41.1566 


56.7201 


5.6682 


44.0711 


0.9796 


0.9938 


0.9957 


10.0304 


1.5 


20 


10 


9.9867 


6.5956 


13.0066 


6.9918 


62.5063 


0.9433 


0.9801 


0.9890 


10.0148 


1.5 


60 


10 


9.3295 


11.0001 


4.9916 


10.5193 


22.7730 


-0.7781 


-0.3566 


0.8522 


10.9880 


0.5 


20 


1.4 


9.9997 


3.3705 


7.8889 


0.7872 


10.6441 


0.9294 


0.9760 


0.9847 


10.0003 


0.5 


50 


1.4 


9.9907 


44.9413 


61.9743 


6.1869 


48.1536 


0.9829 


0.9948 


0.9964 


10.0101 


0.5 


20 


10 


9.9955 


7.2160 


14.2620 


7.6482 


68.5521 


0.9526 


0.9834 


0.9908 


10.0050 


0.5 


60 


10 


9.7394 


12.6647 


8.1550 


11.5177 


32.7242 


-0.5226 


0.0717 


0.8049 


10.3398 


Advanced LIGO 


1.5 


25 


1.4 


9.9995 


1.5726 


2.3120 


0.0769 


2.0304 


0.8002 


0.9397 


0.9319 


10.0006 


1.5 


100 


1.4 


9.6291 


20.2342 


11.4112 


0.3305 


3.9132 


0.9409 


0.9824 


0.9855 


10.4363 


1.5 


200 


1.4 


8.1761 


100.5296 


34.8762 


0.8944 


7.1141 


0.9524 


0.9855 


0.9898 


12.5470 


1.5 


100 


10 


7.6653 


7.6876 


3.6663 


0.8168 


7.2850 


0.4943 


0.8220 


0.8837 


13.5008 


1.5 


200 


10 


4.7707 


16.9594 


2.8159 


1.2761 


3.8796 


-0.6381 


-0.0433 


0.7703 


23.1557 


1.5 


100 


50 


4.1663 


6.7824 


0.9771 


2.2953 


4.5952 


-0.8306 


-0.6018 


0.8899 


28.3348 


1.5 


175 


50 


4.8136 


16.9720 


1.2004 


5.3227 


8.9546 


-0.9283 


-0.9057 


0.9909 


47.1507 


0.5 


25 


1.4 


9.9998 


1.5726 


2.3121 


0.0769 


2.0304 


0.8002 


0.9397 


0.9319 


10.0002 


0.5 


100 


1.4 


9.8716 


20.4346 


11.5270 


0.3338 


3.9530 


0.9421 


0.9827 


0.9858 


10.1475 


0.5 


200 


1.4 


9.2305 


132.9684 


46.3067 


1.1798 


9.4457 


0.9727 


0.9917 


0.9942 


10.9152 


0.5 


100 


10 


8.9967 


11.4032 


5.8388 


1.1846 


11.6961 


0.7646 


0.9230 


0.9461 


11.2882 


0.5 


200 


10 


6.8066 


19.8847 


4.4807 


1.4050 


5.9807 


-0.2786 


0.3719 


0.7714 


15.6651 


0.5 


100 


50 


5.9767 


7.0260 


1.2240 


2.3409 


5.8647 


-0.7460 


-0.3025 


0.7993 


18.2835 


0.5 


175 


50 


5.4400 


17.2493 


1.3221 


5.3654 


9.1950 


-0.9194 


-0.8559 


0.9743 


28.4206 


VIRGO 


1.5 


20 


1.4 


9.9894 


1.5871 


2.8676 


0.1339 


3.0652 


0.8619 


0.9498 


0.9679 


10.0116 


1.5 


100 


1.4 


9.8253 


71.2674 


49.3691 


2.4945 


19.7938 


0.9747 


0.9923 


0.9946 


10.1916 


1.5 


20 


10 


9.8943 


2.7362 


4.0693 


1.0874 


15.4957 


0.8870 


0.9570 


0.9774 


10.1154 


1.5 


50 


10 


9.6406 


6.5355 


6.0863 


2.1034 


19.5048 


0.7740 


0.9210 


0.9537 


10.4273 


1.5 


100 


10 


8.9350 


11.7992 


5.2447 


3.0553 


12.5793 


-0.2535 


0.3875 


0.7797 


11.4628 


1.5 


70 


50 


8.0194 


10.1883 


1.6360 


9.0353 


15.5288 


-0.9005 


-0.7771 


0.9528 


15.1284 


1.5 


90 


50 


8.8495 


12.9812 


1.4506 


11.9826 


19.7503 


-0.9062 


-0.8606 


0.9839 


18.3197 


0.5 


20 


1.4 


9.9964 


1.5874 


2.8682 


0.1339 


3.0658 


0.8619 


0.9499 


0.9679 


10.0039 


0.5 


100 


1.4 


9.9386 


81.8448 


56.7624 


2.8626 


22.7583 


0.9808 


0.9942 


0.9959 


10.0643 


0.5 


20 


10 


9.9641 


2.8685 


4.2790 


1.1410 


16.3023 


0.8973 


0.9610 


0.9795 


10.0386 


0.5 


50 


10 


9.8721 


9.4181 


9.0312 


3.0223 


29.0276 


0.8908 


0.9627 


0.9778 


10.1444 


0.5 


100 


10 


9.6003 


15.8131 


8.7251 


3.8272 


20.7374 


0.2455 


0.7029 


0.8553 


10.5102 


0.5 


70 


50 


8.8855 


11.2685 


2.2537 


9.7760 


18.2478 


-0.8595 


-0.5971 


0.9003 


11.9564 


0.5 


90 


50 


9.0797 


15.2788 


1.9791 


13.9267 


22.9881 


-0.9111 


-0.8146 


0.9655 


13.3618 
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Appendix A. The waveform derivatives 

In order to calculate the Fisher matrix with respect to the {In^iNs, /o^c (f'c, InA^, 
Inr?, Ina, 70, to} basis, we need to compute first the waveform derivatives. For the 
inspiral waveform /iins(/) these are the following 



/ilNS (A.l) 

I 2^(///o)/i:ns (A.2) 

- i hms (A. 3) 

-Y^Ak ik~5)u''-'h,.s (A.4) 



d In Aws 

dhms 

dfotc 
d(j)c 

d\nM " 128 

^ -i— ySfe^^'^-'Vs, (A.5) 
91n?7 128^ /.INS, y ) 

where the parameter Ak are given by equations H25() - (|29|l and Bk are 

Bo = p^=0 (A.6) 
a In 77 

Si = ^ = (A.7) 
o In 77 

dA2 , 743 11 , _2/5 

^^^ah^^(-378 + y^)'^ ' (^-'^ 

omr] 5 

dAi . 3058673 5429 617 2^ -4/5 

Ba = - — ^ = ( 1 77 H r]^)n ^'^ . (A. 10) 

dim] ^ 127008 504 ' 12 ' ^ ' ^ ' 

The inspiral waveform has no dependency on Ina, 70 and to- Therefore the remaining 
derivatives vanish 

d In a djo dto 

The derivatives of the ringdown waveform /ird(/) can be computed by taking into 
account the implicit dependencies of AnuiAms, M, rj, a), /qnr(A^, '7, a), and Q{a). We 
get 



91n Ans 



/iRD (A.12) 



dKn _1~ dKn /A T^^ 

ainX "6'^'^° 91n/Q«R ^"^'^^^ 
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dlnrj 
9 In a 



where 



2 ^ ^ 3 i9feRD 
5 5 9 In /qnr 

9(-100 + 21(l-a)3/i0a) ~ 
40(-100 + 63(1 - a)3/io)(-l + a) ^ 

189a a/ino 
10(-63 + 100(1 - a)Vio + 63a) 9 In /q^r 

9a 

20(1 - a) dlnQ 

^rdc jqnr 



(A.14) 



(A.15) 



91n/QNR 



27r 

e''^«(2iQ + l)Q 



(A.16) 



e-'^"(l -2iQ)Q 



. (/QNR(i - 2g) + 2/Q)2 (2/Q + /q„r(2Q + i)f 

A 2if7rto f 



d\nQ 



2tiQ 



I 



-*7o 



(A.17) 



\(/qnr (2i + ^) -2i/)^ 
The remaining derivatives are 



2^ 



(A.18) 



.*7o 



-170 



^'q'^ " 2i(/ /qnr) '^'q"^ ^ 2i(/ + /qnr)^ 



9io 

9/ird 9/ird 



2i/7r/iRD 
= 



(A.19) 
(A.20) 



Appendix B. Analytical analysis of the Fisher matrix for the ringdown 
waveform 

In what follows we are interested in finding an analytical approximation to the Fisher 
matrix for the ringdown waveform in order to compare and verify the numerical results 
obtained with a Fortran code. For this comparison let us focus with the simpler 
case with five parameters, in which we are interested in computing the Fisher matrix 
with respect to the basis (In^iNSj fotc, (/)c, In-A^, In??), thus assuming that there are no 
uncertainties in the parameters (In a, 70, io)- Of course, the ringdown signal does not 
depend on fotc and (j)c, therefore the problem is reduced to three parameters, although 
the signal would depend only on two independent ones, e.g., (In ^rd, In /qnr). 

The way we proceed is to compute first the Fisher matrix of the ringdown signal 
with respect to (Iu^rd, In /q^r). Assuming constant noise over the bandwidth of the 
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signaljj, and taking the (5-function approximation as in |35| . the elements Ti^a in a ^-J^id 
Tin /qnr In /qnr '^^^ be computcd analyticaUy using Mathematica. For 70 = we get^: 

2Al^Q{l + 2Q-') 

In^R. In^RO - ^ 4Q2)/^^^5„(/qnr) ' ^ 

and 

_ ^LQ(i + 4g^ + 8g«) 

in /QNR in /QNR ^ 4Q2)/^^^5„(/^^^) ■ ^^-^^ 

For the cross term ri,i_4 in /qnr > Finn's approximation ^5] can no longer be employed, 
because dh^o/dln /qj^^ is not a symmetric function around /qnr- In order to compute 
this term we will consider the following properties we have derived. 

The reader should notice that for any set of parameters (In^, {A*}) and any 
waveform of the form 

h{AAXlJ)^AH{{y}J), (B.3) 

the elements of the Fisher matrix satisfy the relations 

hnAX' ~ 2 ' ^ ' 

T^J^^J^lnAX'-{h\^^Jh). (B.5) 

Then using the standard definition of SNR given by Eq. ^ we have 

^InAinA^P^, (B.6) 

and consequently 

Tln^A. = (B.7) 

These relations hold true for both the inspiral and the ringdown signals when 
considering A to be the amplitude of the signal. In case of the ringdown signal, 
using equations (|B.ip . ljB.6p and (|B.7p we get 



2Al^Q{l + 2Q^ 

NR ) 



Pro - , f^2\Y- ETTT T ' (^•' 



where 



Tln^RD In /qnr - ^2^""° (1 + '5') , (B.9) 

S= .,l dS(f ) ^ (B.IO) 
'-'jiv/qnrJ " in Jqnr 

Note that the (5-function approximation in this case is equivalent to consider 5 = 0, 
but this term is not negligible. For example if we consider initial LIGO and a total 
mass of 10, 20 or 100 Mq, the corresponding S value would be 1.989, 1.959 and 1.264 
respectively. 

II The approximation that the noise is constant over the bandwidth of the signal is a good 
approximation for all the detectors considered here when a > 0.9 corresponding to A///qnr < 0.5 
as explained in 1351 . In this paper we consider only the case in which a = 0.98. 

^ If instead of using 70 = wo take 70 = 7r/2 then rin_4j^j^ inylao becomes 4^^^ Q''/[7r(l + 
4Q'^)/qnrSti(/qnr)] equivalent to Finn's result 
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The Fisher matrix with respect to the basis (In A^s , hi , hi rj) (which naturaUy 
wiU be degenerate) can be easily be computed by taking into account the amphtude 
relation given by Eq. H45|) and considering 

dhu.u/d\nAms \ / 1 \ / /I] A \ 

dK,^/dlnM =1/6 -1 7° .(B.ll) 

dK^/d\nv J V 2/5 3/5 / V Ohr^u/dln}^^^ J 

If we define the constant matrix 

/I \ 

C= 1/6 -1 , (B.12) 
V 2/5 3/5 ; 

let r be the Fisher matrix with respect to (hi^ao, In /qnr) and F the Fisher matrix 
with respect to (In^i^s, In A^, In 77), in this particular case, F and F are related in the 
following way 

f=CFC^, (B.13) 
where the superscript T indicates transposed matrix. The matrix F has the elements: 

Tin^iiNs in^iNs = pId ' (B-14) 
f 2 4 + 35 

i^lnAmslnM - Pro q ' (aAb) 

TinANsii") = Prd — — ' (B.16) 
2 2 (72Q2 + 65 + 43) g2 + 65 + 25 

T,nMlnM=PRu 72Q2 + gg ' (^.17) 

fi„A,i„„ = /9L^ (-72Q2 + 95-^^i^ + 13) , (B.18) 

fin „ In, = ^ (36g' - 125 + - 4^ , (B.19) 

and trivially 

f/ot.A. =f^.A. =0. (B.20) 

The analytical approximation and the numerical results are compared in tablc lBTI 
for the initial LIGO detector and two different values of the total mass. 
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